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Abstract 

Results on J/ip production in ep interactions in the HI experiment at HERA are presented. 
The J /if) mesons are produced by almost real photons (Q 2 « 0) and detected via their 
leptonic decays. The data have been taken in 1994 and correspond to an integrated 
luminosity of 2.7 pb -1 . The jp cross section for elastic J /if) production is observed to 
increase strongly with the center of mass energy. The cross section for diffractive J /if) 
production with proton dissociation is found to be of similar magnitude as the elastic cross 
section. Distributions of transverse momentum and decay angle are studied and found to 
be in accord with a diffractive production mechanism. For inelastic J /if) production the 
total 7p cross section, the distribution of transverse momenta, and the elasticity of the 
J /if) are compared to NLO QCD calculations in a colour singlet model and agreement is 
found. Diffractive if)' production has been observed and a first estimate of the ratio to 
J /if) production in the HERA energy regime is given. 
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1 Introduction 



The high center of mass (cms) energy of \fs ~ 300 GeV available at the electron proton collider 
HERA has led to a renewed interest in the study of J / if) production in lepton proton scattering. 
The production cross section is dominated by photoproduction, i.e. by the interaction of almost 
real photons, Q 2 ~ 0, where Q 2 is the negative four momentum transfer squared to the scattered 
lepton. The topic of this paper is an analysis of photoproduction of J / if) mesons in elastic and 
inelastic processes. 

Several mechanisms have been suggested to describe photoproduction of J /if) mesons 

•yp — > J /if) X. 

For the description of the elastic (or exclusive) process where X = p, a diffractive mechanism 
has been proposed by Donnachie and Landshoff |T[]. The mechanism is based on Regge phe- 
nomenology where J /if) production is mediated by pomeron exchange (Fig. 1(a)) and can be 
extended to include also proton dissociation (Fig. 1(b)). Measurements at HERA are expected 
to shed new light on the diffractive production mechanism and the nature of the pomeron ||. 
At lower center of mass energies pomeron exchange was successfully used in the framework of 
the vector dominance model (VDM) 0, f|] to describe the production of light vector mesons. 
The predicted slow variation with energy of the elastic cross section was confirmed and an 
exponentially falling t distribution was found, where t is the squared four momentum transfer 
to the scattered proton. When applied to J / if) meson production at low center of mass energies 
the prediction of the VDM model was found to be more than an order of magnitude above the 
data but other features like energy dependence and t distributions agreed within errors. 
Attempts have been made to describe elastic J /if) production in perturbative QCD. In the 
approach by Ryskin |5|] the interaction between the proton and the charm quark is mediated 
by the exchange of a gluon ladder (Fig. 1(c)). Non perturbative effects are included in the 
gluon distribution of the proton which enters the production cross section quadratically. A 
measurement of elastic J/ if) production could therefore be a sensitive probe of the gluon density 
in the proton. Recently higher order effects have been calculated in this model and a comparison 
to preliminary HERA data |J was carried out 0. The Ryskin model predicts a different energy 
dependence of the elastic J/ if) cross section than the Donnachie-Landshoff approach. The energy 
dependence in the Ryskin model is coupled to the low x behaviour of the gluon density in the 
proton. Using a gluon density increasing towards low x which describes recent measurements 
of the structure function F2 at HERA||, results in a fast increase of the cross section for elastic 
photoproduction of J /if) 0. 

In elastic J /if) production at small momentum transfer the J /if) meson retains approximately 
the full photon energy (z ~ 1 with z = E^/E^ in the proton rest system). J /if) production 
with proton dissociation although strictly speaking an inelastic process also leads to z values 
close to 1. In contrast, inelastic processes will have z values below 1 and a high mass hadronic 
state is formed. The principal inelastic process proposed is the photon gluon fusion mechanism 
where the photon emitted by the incoming lepton interacts with a gluon from the proton via 
the charm quark (Fig. Id). This mechanism is calculable in perturbative QCD due to the 
hard scale given by the mass of the charm quark. Attempts have been made to determine the 
gluon density in the proton from this reaction in several fixed target experiments |T0|, [TT 



In the colour singlet model [12| for photoproduction of J /if), which is based on this picture, 



the formation of a J /if) state is accompanied by the emission of a hard gluon. Comparing 
predictions of the colour singlet model to data a discrepancy in absolute magnitude was found 
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(c) (d) 

Figure 1: J/ip production mechanisms: (a) elastic J/ip production via pomeron exchange; 
(b) diffractive proton dissociation; (c) elastic J / ip production in perturbative QCD: exchange 
of a gluon ladder ||/; (d) photon gluon fusion model for inelastic J / ip production (colour singlet 
model 0). 



which was attributed to missing higher order calculations ( "K-factor" ) . Subsequently several 
improvements have been proposed which led to better agreement with the data O]. Recently, 



complete next-to-leading order (NLO) calculations have been performed [14, 15] and compared 
successfully to fixed target data and to preliminary HERA data. 

Motivated by TEVATRON data on inelastic J ftp and ip' production [TfJ] which may require 
colour octet contributions to explain the measured cross sections, calculations have been per- 
formed for possible colour octet contributions in photoproduction at HERA [18]]. Both 
groups attempt to estimate the colour octet contributions to J/ip mesons at high z; in [|17] the 
contribution to inelastic J /ip production is estimated. 

We present here an analysis of elastic and inelastic J/ip production in positron proton collisions 
near Q 2 = for 7p cms energies up to 150 GeV. The analysis is based on J/ip decays to leptons, 
J/ip — > n + or -> e + e~. The data were collected with the HI detector at HERA and 
correspond to an integrated luminosity of ~ 2.7pb~\ 

The inelastic process is analysed for the first time in HI; preliminary results have also been 
shown by the ZEUS collaboration H. The analysis of the elastic process is an update of 



a previous letter |T9| where we have presented a measurement of cr(7p — > J/ip + X) which 
showed a strong increase of the cross section with W 1P , the photon proton center of mass 
energy, compared to experiments at lower cms energy. This increase was faster than expected 
from the Donnachie-Landshoff prediction. At that time contributions from processes with 
proton dissociation could however not be excluded completely. With the increased statistics 
now available and an improved analysis method these inelastic processes can be efficiently 
recognised. The fast increase of the elastic cross section with W 1P was also observed by the 
ZEUS collaboration 
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The paper is organised as follows. After a brief introduction of the kinematics (section 2), the 
experimental conditions and the pre-selection of the data sample are described in section 3. The 
analysis of the elastic J/ip samples follows in section 4, which includes the energy dependence 
of the total 7p cross section and the distribution of pf, the transverse momentum of the J/ip. 
The second part of section 4 contains the analysis of the proton dissociation process and the 
decay angle distribution of diffractively produced J / ip mesons. In section 5 the selection of the 
inelastic events is discussed and results are presented including the energy dependence of the 
7P cross section, the distributions of pi and of the elasticity z of the J/ip. 
A first cross section for ip' production in the HERA energy range is given in the last section. 



2 Kinematics 



The variables used for the description of J/ip photoproduction at HERA are: 



Q 2 = ~q 2 = -(k-k'f (1) 
t = (P p -Pf (2) 
s = (P p + k) 2 (3) 

(4) 



P P -q 



Pp ■ k 



W 1P = ^(P P + q) 2 (5) 

with the 4-momenta P p and k of the incoming proton and positron, k! of the scattered positron, 
q = k — k! of the exchanged photon, and P' of the system X, which is identical to the proton 
for elastic J/ip production. 

The Bjorken variable y can for photoproduction be approximated by y ~ E^/E e , where E 1 and 
E e are the energies of the exchanged photon and the incoming positron. The variable y can be 



computed from the observed final state using the method by Jacquet-Blondel [21 



= Z(E- P z ) = (E-P z ) J/lP + Z iest (E-P z ) 

2E e 2E e ' U 

The sum is over all visible particles, i.e. tracks and energy deposits in the calorimeter. In 
order to avoid double counting the calorimeter energy in a cylinder of radius 30 cm around the 
extrapolated charged tracks is excluded. For events where only the two J/ip decay leptons are 
observed in the tracking detector the calorimeter information is not used for the calculation of 
y. Analysing photoproduction it is customary to use W 7P : 



7 2 p = ys-Q 2 + m 2 p 



The separation of the events into an elastic and an inelastic sample utilizes the topology of the 
event and the variable z which is defined as: 

P P 

= ^jl_L± (7) 

Pp ■ q K ) 

where P^ is the four vector of the J/ip. Using eqn. (4) z can be expressed as: 
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- with _ , (8) 



For elastic events y^ = y and thus z = 1. 

The ep cross section and the 7p cross section are related by: 

a(ep^eJ/^X)= dy dQ 2 ■ f j/e (y,Q 2 ) ■ a rp (Q 2 ,y) (9) 

(y) 

where the cross section on the left side is the measured ep cross section covering a range 
Qmin — Q 2 — Qmax an d determined in bins of y. The kinematical minimum for a given value 
of y is Qmm, and Q 2 max is the effective upper limit of the selected event samples: 



Qmin ~ m e -i „ ! Qmax ^ GeV 



The flux of transverse photons is [22 



The longitudinal photon flux amounts to 2% of the transverse flux for the given kinematical 
range. Since the photon flux decreases rapidly with Q 2 and y, the weak Q 2 and y dependences 
of a 7 * p (Q 2 ,y) play only a minor role and the photoproduction cross section ai^p — > J/ipX) 
can in first approximation be identified with cr 7 * p (Q 2 , y). It is then obtained as: 

a(jp -»■ J/tpX) = a(ep ^ e J/tp X)/% /e 

where $ 7 / e is the photon flux integrated over Q 2 and y. The correction for this approximation 
is small; it can be taken into account as a small change of the W w value at which the cross 
section measurement is performed. 



3 Experimental Conditions 

The data were taken in 1994 with the HI detector operating at the electron proton storage 
ring HERA, where positrons of 27.5 GeV collide with protons of 820 GeV. In 1994 HERA was 
operated with 153 colliding positron and proton bunches. The integrated luminosity used for 
this analysis is 2.7pb _1 for the decay J/tjj — > and 2.0 pb -1 for J/ip — > e + e~. The HI 

detector is described in P5| . We repeat here the essential features of the detector components 
used for the analysis. 
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3.1 The HI Detector 



The central tracking system is mounted concentrically around the beamline and covers polar 
anglesQ between 20° and 160°. Measurements of charge and momenta of charged particles are 
provided by two coaxial cylindrical drift chambers (central jet chambers, CJC) ||24j| . Two sets 
of cylindrical drift chambers for measurement of the z-coordinate and multiwire proportional 
chambers (MWPC) for triggering are placed at two radial positions. One set surrounds the 
beamline within the inner CJC and the other is mounted in between the two jet chambers. 
The central tracking system is complemented by a forward tracking system which covers polar 
angles 7° 55 6 ^ 25°. In the present analysis the forward tracker is only used to detect events with 
tracks other than the J/ip decay leptons. 

The tracking system is surrounded by a highly segmented liquid argon (LAr) sampling calorime- 
ter pS| with an inner electromagnetic section consisting of lead absorber plates with a total 
depth of 20 to 30 radiation lengths and an outer hadronic section with steel absorber plates. 
Polar angles between 4° and 153° are covered by the calorimeter and the total depth is 4.5 to 8 
interaction lengths depending on the polar angle. The backward region, 155° ^6*^176°, is cov- 
ered by a lead scintillator calorimeter where the scattered positron is detected for Q 2 ^4 GeV 2 . 
The magnetic field of 1.15 T is produced by a superconducting solenoid surrounding the LAr 
calorimeter. 

The iron flux return yoke surrounding the superconducting solenoid is instrumented with lim- 
ited streamer tubes to provide muon identification; it is segmented into 10 iron plates of 7.5 cm 
thickness and instrumented with up to 16 layers of streamer tubes. Muon tracks are recon- 
structed in the region 4° <9 < 171° with a spatial resolution of the order of 1 cm. 
In addition to this central muon detector there is a toroidal muon spectrometer outside the 
main HI magnet covering small polar angles. In the present analysis its driftchambers which 
cover 3° 17°, are used for the recognition of events with proton dissociation. For the same 
purpose a system of scintillators - the proton tagger - is placed 24 m downstream the proton 
beam around the beampipe, covering an angular range of approximately 0.06° ^ 9 ^ 0.25°. 
The luminosity is measured using the radiative process ep — *> ep~y where the photon is detected 



in a luminosity monitor [26, 23 



3.2 Trigger and Data Processing 

The background rate at HERA is high, mainly due to interactions of the beam protons with 
gas in the beam pipe or surrounding material and due to photoproduction of light quarks. 
Therefore a restrictive trigger is necessary which cannot use the distinctive signature of the 
scattered positron since for photoproduction the positron stays in the beampipe^]. The trigger 
essentially has to rely on the decay leptons of the J/ip. 

Compared to our previous report |19| the trigger for muons and for electrons from J/ip decays 
has been improved resulting in an approximate efficiency for electron pairs of 50% and for muon 
pairs of nearly 60% with a tolerable background rate. 
The following triggers are utilized: 

■"■HI uses a right-handed coordinate system defined as follows: the origin is at the nominal interaction point 
with the z-axis pointing in the proton beam direction, hence the polar angle 9 is measured with respect to the 
proton beam direction. The region of small polar angles is called "forward" . The plane perpendicular to the z 
axis is named r — <f> plane. 

2 The scattered positron is with a few exceptions not detected in the low angle tagger, because there was 
practically no overlap of its sensitive region (y £ 0.3) with the bulk of the J/tp data which are at y $ 0.25. 
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1. Track triggers from Multi Wire Proportional Chambers (MWPC) |27J demand- 
ing the origin of the tracks in the z direction to be near the nominal vertex region. 

2. A coplanar track trigger constructed from the MWPCs which demands exactly two 
tracks of transverse momentum p t Z 500 MeV, roughly coplanar with the beams. 

3. A driftchamber track trigger demanding one track in the CJC with a trans- 
verse momentum p t > 450 MeV which originates from the interaction point in the plane 
perpendicular to the beams within ~ 2 cm. 

4. Muon triggers demanding a penetrating particle detected in the central muon systems. 

5. A calorimetric trigger for low energy electromagnetic clusters with E > 800 MeV 
roughly aligned with a track candidate in the proportional chambers. 

These triggers are combined such that each event class is triggered by at least two different 
trigger combinations which can be compared with each other for determination of the efficiency 

For exclusive J/ip meson production in both leptonic decay modes a combination of elements 
1, 2 and 3 gives a trigger based purely on tracking chambers. A trigger for muons - which is 
used for elastic and inelastic processes - is based on track triggers 1 and 3 and the muon trigger 
4. For electrons the combination of track triggers 1 and 3 with the calorimetric trigger 5 is 
used imposing a requirement of low track multiplicity in the MWPCs therefore only sensitive 
to elastic and proton dissociation processes. 

3.3 Track Selection and Lepton Identification 

The event selection starts from tracks found in the central drift chambers CJC which have 
been associated to the primary e + p interaction point (vertex) in a constrained fit which helps 
to increase momentum resolution and to reduce pattern recognition ambiguities. 

Electrons are identified in the electromagnetic section of the LAr Calorimeter by linking a 
reconstructed drift chamber track to a calorimeter cluster with energy E duster > 0.8 GeV and 
demanding the measured energy to be comparable with the momentum [Eduster I Vtrack > 0.7). 
The efficiency for electron identification is ~ 85% for p > 0.8 GeV. 

A particle is identified as a muon if the track in the drift chamber is either linked to a track 
element reconstructed in the central muon detector or if it is identified as a minimum ionizing 
particle in the LAr Calorimeter. For the link with a track in the central muon system a 
drift chamber track is extrapolated taking into account the deflection in the magnetic field, 
the energy loss in the material of the detector and multiple scattering. The \ 2 probability 
comparing the parameters of the two tracks is required to be above 1%. A muon signature 
in the LAr Calorimeter is defined by an energy deposit below 3 GeV around the extrapolated 
track, at least three active cells in the hadronic part of the calorimeter, and the particle has 
to penetrate at least 90 % of the calorimeter. With these requirements the thresholds are 0.8 
GeV for the identification of muons in the LAr calorimeter (« 75% efficiency above ~ 1 GeV) 
and 1.5 GeV in the central muon system (« 80% efficiency above ~ 2 GeV). 
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3.4 Event Simulation and Efficiency Determination 



Monte Carlo (MC) techniques are used in order to determine the geometrical acceptance, the 
trigger and selection efficiencies. Events are generated using models and the detector response 
is simulated in detail. The simulated events are subjected to the same reconstruction and 
analysis chain as the data. 

Two models are used to simulate J / ip production. Inelastic J/ ip production is simulated by the 



generator EPJPSI |29j which is based on the colour singlet model in leading order. Diffractive 
J/ip production - elastic and proton dissociation - is modelled by the Monte Carlo generator 
DIFFVM [0. 

Diffractive events are generated with an energy dependence proportional to W^ p 9 , an exponential 
t distribution ~ exp (— b\t\) with a fixed slope of b = 4 GeV -2 in the elastic mode and b = 
2 GeV -2 for proton dissociation. The proton dissociation mode has an additional parameter, 
the mass Mx of the dissociated proton. Events are generated with a \/M\ distribution for 
M x above 4 GeV 2 . Below 4 GeV 2 the distribution is closely modelled to diffractive proton 
dissociation data obtained from deuterium measurements |JT|. The system X is fragmented at 
masses above 2 GeV by treating it as a system of quark and diquark and using the Lund string 
model El. Below 2 GeV it is treated as a nucleon resonance and decays into nucleon and pions 
according to the known branching ratios. 

The detector response in the Monte Carlo simulation is checked in detail by comparing to the 
data. After applying small overall correction factors the dependence of the efficiencies for track 
reconstruction, particle identification and triggering on the polar angle and the momentum of 
the J/ip decay particles is well described by the simulation. This tuned Monte Carlo simulation 
is then used for correcting the data and remaining differences between data and Monte Carlo 
are included in the systematic errors. 



3.5 Pre— selection of the J/ip Samples 

Elastic and inelastic J/ip events in general look quite different in the detector, elastic events 
having only the decay leptons measured while inelastic events are characterized by additional 
tracks. Initially a common selection aims to identify lepton pairs in the J/ip mass region 
irrespective of any other activity in the detector. Two tracks are selected with a common origin 
at the beam interaction point in the r — cp plane and momenta above 0.8 GeV in the range of 
polar angles 20° < 6 < 160°. Both tracks have to be identified as muons or electrons. The 
^-coordinate of the event vertex has to be within 40 cm of the average beam collision point. 
For low multiplicity events, one of the main backgrounds is due to cosmic ray muons, which 
are efficiently rejected by demanding the angle between the lepton candidates to be less than 
177°. 

Photoproduction events are selected by requiring no scattered positron be visible in the 
calorimeter, i.e. below a polar angle of 176°. Rejecting energy clusters above 8 GeV restricts 
the photon virtuality to Q 2 ;$4 GeV 2 . 



4 Analysis of the Elastic and Proton Dissociation Pro- 
cesses 

Starting from the pre-selected data a two track sample is selected for the determination of the 
cross section for elastic J/ip production and for proton dissociation processes. The selected 
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Figure 2: Mass distribution for u. + u~ (a) and e + e~ pairs (Tjj for the two track selection above 
2 GeV. The curves are fits of a Gaussian plus a polynomial background to the J j ip mass 
region. The shaded histogram shows the contribution of QED lepton pairs. For fi + fi~ the 
maximum of the fit is at (3.10 ± 0.01) GeV with a width of 76 MeV. For e + e~ the maximum is 
at (3.08 ± 0.02) GeV and the width is 77 MeV. For both cases the detector simulation yields a 
width of 65 MeV. 



lepton pairs are required to be the only tracks coming from the beam interaction point in the 
sensitive region of the tracking detectors (7° ^8^,165°). 

The reconstructed invariant mass of the lepton pairs for these events is shown in Fig. |2|a) and b) 
for electrons and muons, respectively. A cut around the nominal mass \m^+^ — m^\ < 225 MeV 
is applied. This yields about 400 (350) J/ip candidates for muon (electron) pairs. 

As can be seen from Fig. ^ the background below the J/ ip mass peak for this two track sample 
is low. It can be described by lepton pairs produced via a two photon process, where a photon 
is emitted by the incoming positron and by the proton. This process is calculable in QED. The 



calculation is implemented in the Monte Carlo generator LPAIR |j33| , the result of which is also 
shown in Fig. |2|. 

A subtraction is made for the background below the resonance, approximately 5% for J/ip — > 
and 12% for J/ip — ► e + e~ . The number for electron pairs is larger because it contains two 
effects in addition to the QED lepton pairs. There is an additional background of misidentified 
hadrons with a decreasing mass spectrum which has to be subtracted; on the other hand there 
is a small loss of signal events due to the low mass tail of the J/ip peak. The latter is due to 
radiation in the material of the detector according to the detector simulation. 

At this stage the two track event sample is composed of elastic events and those with proton 
dissociation where the fragments of the excited proton are not visible in the tracking detector. 
The events are then classified into two samples, one with and one without activity in the forward 
detectors. The sample with forward activity will be used in section |4.2j to determine the cross 
section for J/ijj production with proton dissociation. For elastic J/ip production the sample 
without signals in the forward detectors is used. 
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The forward detectors close to the proton beam (see section 34.) are the proton tagger, the 
drift chambers of the forward muon spectrometer, and the low angle region of the liquid argon 
calorimeter (4° & 9 $ 10°). They are sensitive to particles originating from secondary interactions 
of the excited nucleon fragments with material of collimators and the beampipe. Their efficiency 
has been studied with data p4f . 

The three detector systems have different thresholds as regards the mass Mx of the dissociated 
nucleon system X . The proton tagger is sensitive down to the lightest nucleon excitation, 
the forward muon detector is sensitive for 1.5 GeV and the low angle region of the LAr 

calorimeter for Mx > 3 GeV. 



4.1 Elastic J/ift Production 

Monte Carlo studies show that after the two track selection, 85% of the events with proton 
dissociation are recognised in the forward detectors. The elastic event sample is selected by 
rejecting those events containing an energy deposit above 0.75 GeV in the forward LAr calorime- 
ter, or having hits in the proton tagger or having more than one hit pair in the forward muon 
system. The event numbers are listed in table [l]. 

The remainder of the proton dissociation background is subtracted statistically using the Monte 
Carlo simulation program DIFFVM of J/ip production with proton dissociation. The simulated 
events are normalized to the number of events with forward detector signal. All three forward 
detectors independently lead to the same normalization within errors. In table |] the correction 
factor which is applied to account for remaining proton dissociation background is indicated in 
the line marked T — fp.diss? ■ 

Additional background from if}' decays is removed taking into account the measured fraction 
of ip'/ijj production of ~ 20% (see section an d applying the selection cuts to the simulated 
cascade decays into J /ip+anything. A correction of 4 ± 2% is applied. 

The remaining event sample has been examined by a visual scan in order to detect unrecognised 
background. A small number of events has energy deposits in the calorimeter which might 
indicate a background. This number is less than the background expected due to ip' decays 
and no further correction is applied. 



Cross section as function of W 7P 

The elastic ep cross section is evaluated in bins of W lp in the range of 30 - 150 GeV. The Q 2 
range is limited to Q 2 ^4 GeV 2 with an average of 0.13 GeV 2 . The cross section was calculated 
according to: 

Nbg.corr. (1 f p. diss.) (1 ftp') 
O'ep = ~ 
^sel ^trigger ^acc 

where the factors are the number of events corrected for background below the mass peak 
Nbg.corr.) the background fractions of proton dissociation and tp' f p .diss. and fy, respectively, the 
efficiencies for data selection e se i (including track reconstruction, selection cuts, and lepton iden- 
tification) efficiency for triggering trigger-, the geometrical acceptance e acc , and the integrated 
luminosity C The results for the ep cross section are given in table [I] taking into account the 
branching ratio for one leptonic decay channel of (6.0 ± 0.25)%. The ep cross section can be 
converted into a jp cross section using the relations of section 0. The integrated photon flux 
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W 7P [GeV] 


30-60 


60-90 


90-120 


120-150 


J /if} -> n + V 


events 

(bg. corrected) 


52.0±7.5 


64.0±8.6 


80.0±10.0 


37.0±6.6 


€(LCC 

1 fp.diss. 

^selection 

^trigger 


0.512±0.020 

0.88±0.12 
0.486±0.039 
0.371±0.030 


0.844±0.034 

0.87±0.11 
0.396±0.032 
0.490±0.039 


0.669±0.027 

0.88±0.11 
0.485±0.039 
0.724± 0.058 


0.390±0.016 
0.87 ±0.11 
0.581±0.046 
0.828±0.066 


a(ep — > eJ/^p)[nb] 


2.9±0.4±0.5 


2.0±0.3±0.4 


1.7±0.2±0.3 


1.0±0.2±0.2 


J/t/j — ► e+e 


events 

(bg. corrected) 


48.4±7.5 


62.5±8.5 


43.1±7.1 


10.6±3.5 


^seZeciion 
^trigger 


0.550±0.038 
0.474±0.033 


0.515±0.036 
0.550±0.039 


0.447± 0.031 
0.684± 0.048 


0.223±0.016 
0.791±0.055 


<r(ep — > eJ/^p)[nb] 


2.5±0.4±0.4 


1.8±0.2±0.3 


1.5±0.2±0.3 


1.1±0.3±0.2 



Table 1: Elastic J/ij) production: Efficiencies and cross sections for elastic J/ip production 
measured via the leptonic decays to fi + fi~ and e + e~ are given in intervals of W 1P (Q 2 5= 4 GeV 2 ) 
in the acceptance region. The number of events (two track selection requiring no forward 
tag) is corrected for background below the mass peak. The remaining proton dissociation 
background fp.diss. is calculated by Monte Carlo. The acceptance e acc describes the probability 
for both leptons to be in the 8 range 20° < 6 < 160°. The selection efficiency e selection contains 
the efficiencies for track reconstruction, lepton identification and selection cuts. The errors of 
acceptance, proton dissociation background, selection and trigger efficiency give the systematic 
errors. The first error of the ep cross section is statistical, the second one systematic. 



j /if) -> m- 


W lp [GeV] 


30-60 


60-90 


90-120 


120-150 


W [QcV] 


42 


72 


102 


132 




0.0736 


0.0370 


0.0229 


0.0149 


v(lP^ J/#) Wi [nb] 


39.6±5.7±7.2 


53.4±7.2±9.8 


76.2±9.5±13.9 


67.3±11.9±12.3 


a(^p — > J / ipp) ee [nb] 


34.4±5.3±6.0 


48.2±6.6±8.4 


63.8±10.4±11.1 


71.0±23.4±12.4 


cr(7p — * J/ipp)[nb] 


36.8±3.9±6.6 


50.6±4.8±9.1 


70.6±7.0±12.7 


68.0±10.6±12.2 


<j{lP^ J/ipX)[nb] 


23.0±3.2±4.0 


63.5±5.8±11.4 


62.7±7.4±11.2 


128.9±19.5±23.2 



Table 2: Elastic J/ij) production and J/ip production with proton dissociation: 7p cross 
sections for the elastic production of J / ip mesons are given for decays to muons and electrons 
separately and also combined. In the last line the combined cross section for J/ip production 
with proton dissociation is given. W is the corrected bin center. $ 7 / e is the photonffiix 
integrated over Q 2 and W 1P . The first error of the cross section is the statistical error, the 
second one systematic. 
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is shown in table 0. The center W of the W 7P bins of 30 GeV width is calculated taking into 
account a y dependence of the 7p cross section corresponding to ~ W® p 9 , a Q 2 dependence of 
<7 7 *p as given by the VDM propagator (m 2 /(Q 2 + m 2 )) 2 , and the longitudinal cross section to 
be Q 2 /m 2 times the transverse cross section. The difference to the uncorrected bin center is 
estimated to be —3 GeV. 

Since the data based on J /if) decays into ji + n~ and e + e~ agree within errors the 'jp cross 
sections have been averaged. The resulting combined cross section is shown in Fig. [3] and listed 
in table The data show a rise with W 1P which in the HERA regime can be represented as 
<7 7P oc W^ p . A fit to the HI data yields 5 = 0.64 ± 0.13, where the error includes statistics 
and systematical effects. In Fig. |3| also results by the ZEUS collaboration [20| in the same 
energy range as the present data are shown, which agree well. Using also the ZEUS data and 
the measurements at lower energies [[35|, |36|, |37| which are also shown in Fig. |3|, a combined fit 
yields 5 = 0.90 ± 0.06. 

The energy dependence expected from the Donnachie-Landshoff model |1[ corresponds to 5 = 
0.32 (neglecting any shrinkage, see also next section) and is also shown in Fig. |^. It falls below 
the HERA data by more than a factor of 3 corresponding to more than 3 standard deviations 
if normalised at fixed target energies. 

The prediction of the QCD model due to Ryskin and including higher order corrections |7J 
is also compared with the data in Fig. |3|. The prediction depends quadratically on the gluon 
distribution taken at the scale (Q 2 + m 2 b )/A m 2.4 GeV 2 , which can be parametrised as x g(x) oc 
x~ x at low values of x, the fraction of the proton momentum carried by the gluon. For the 



present data x ~ mLjW 2 ~ 10 3 . Using the gluon distribution function from MRSA' [38 



corresponding to A ~ 0.2, gives good agreement with the data between 10 and 150 GeV. The 
sensitivity to the gluon distribution is illustrated by comparison to the GRV fit |39| which 
corresponds to A ~ 0.3 — 0.4. The corresponding curve in Fig. |3| with the parameters used in |7J 
yields a steeper energy dependence than the data. 



Systematic errors 

A breakdown of the systematic error of the cross sections in tables [l] and |2| is given in table |||. 
The main contribution to the systematic error of elastic and proton dissociation cross sections 
is due to the separation of these two event classes. The uncertainty is estimated to be 12% 
by varying the cuts on the three forward detectors and by using different combinations of two 
of them. The uncertainty of the background which was subtracted statistically in the elastic 
cross section was estimated to be 2% and is included in the total of 12%. It was estimated by 
varying the assumed \ jM\ dependence of the cross section from k=2 to k=2.5 (k = 2.2 ± 0.2 
was measured in fjOfl). 

The systematic errors of the track reconstruction efficiency (3%), single muon identification 
(5%), electron identification (4%) and trigger efficiency (9% and 8% for muon and electron 
triggers) are estimated by comparing the Monte Carlo efficiencies with the efficiencies deter- 
mined from the data and using remaining differences as error. 

The systematic error of the angular acceptance depends mainly on the energy dependence of 
the 7p cross section and is estimated by varying a lp oc between 0.6 < 8 < 1.0 which results 
in a 4% error. The uncertainty in the background due to ip' cascade decays via J/ip accounts 
for the measurement error on the ratio of ip'/ip production. For the evaluation of the 7p cross 
section an uncertainty in the photon flux of 2% is estimated by varying the upper limit of the 
Q 2 integration by ±1 GeV. 
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Figure 3: a) Total cross section for elastic J/ip photoproduction. The inner error bars of the 
HERA points are statistical, the outer ones contain statistical and systematic errors added in 
quadrature. The data at lower cms energies are from previous experiments |3^, |36|, [3/]/; they 
were corrected with the new J / if) decay branching ratio where necessary and include systematic 
errors (added in quadrature). A parametrisation of the energy dependence as W^ p with 5 = 0.90 
(full curve) and 5 = 0.32 (long dashed) normalized to the E 516 point is shown. Also shown is 
the result of calculations according to the Ryskin model FX) using the MRSA ' (dashed) and GRV 
parametrisations (dash dotted) of the gluon density, b) Total cross section for J / -production 
with proton dissociation. The inner error bars of the HI points are statistical, the outer ones 
contain statistical and systematic errors added in quadrature. The full line represents a fit to 
W* p with 5 = 1.2 ±0.2. 
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H + H e + e 


Separation elastic/p.diss. 


12% 


Trigger 


9% 8% 


Single lepton identification 


5% 4% 


Single track reconstruction 


3% 


Acceptance 


4% 


Branching ratio 


4% 


Photon flux 


2% 


Luminosity 


1.5% 


if)' -background 


2% 


Total 


18 % 17% 



Table 3: Systematic errors 



■pi distribution 



The t-distribution is one of the distinctive features of a diffractive process. It is expected to be 
exponential and the slope parameter in Regge type models is expected to increase with energy 
(shrinkage) as b = bo + 2a' In W 2 /Wq with a' = 0.25 GeV~ 2 for pomeron exchange In 
contrast, models based on perturbative QCD predict little shrinkage. 

Since the scattered electron is not measured, t can only be approximately determined for the 
data assuming Q 2 = 0, then t ~ —p\ where p\ is the transverse momentum of the J /if} with 
respect to the beam axis. The data analysis was carried out as in the previous section, but 
the efficiencies and background proportions were determined in bins of p t . The cross section 
da ep /dp 2 is shown in Fig. f| a. A steep slope is observed for pj < 1 GeV 2 with a tail towards 
higher values. The tail is compatible with Monte Carlo expectations for events with Q 2 ^ 0. 

The slope obtained by a log-likelihood-fit of an exponential exp (—bp 2 ) to the data below 
p 2 < 1 GeV 2 is 

b = (4.0 ±0.2 ±0.2) GeV" 2 . 

The first error is statistical. The second one contains the dominant systematic contributions 
estimated by varying the upper limit of the fit region in p\ between 0.75 and 1.25 GeV or by 
using a x 2- fit, alternatively. The error in the slope parameter due to using p 2 t instead of t is 
calculated in the simulation to be —10%. 

For the W yp range from 30 GeV to 90 GeV the fit yields b = (3.7 ± 0.3 ± 0.2) GeV" 2 and 
b = (4.5 ± 0.4 ± 0.3) GeV~ 2 for the interval from 90 GeV to 150 GeV. Within the errors no 
clear evidence for shrinkage is observed. The expected change in the b slope from the Regge 
prediction is only 0.6 GeV -2 in the HERA energy range. This is of the same order of magnitude 
as the experimental errors, therefore no conclusion can be drawn with present statistics from 
the data in the HERA energy range alone. 

The experimental situation at low cms energies is unclear: 6-values ranging from b ~ 3 GeV~ 2 
to b ~ 5 GeV -2 were measured || [K| |36|, fll], (|!| thus also preventing any conclusion about 
shrinkage. 
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Figure 4: a) da /dpi for elastic J/ip production integrated over 30 GeV < W yp < 150 GeV. 
The error bars contain statistical and systematic errors added in quadrature. The straight line 
is a fit to the data in the range p 2 < 1 GeV 2 of the form exp(—bp 2 ) with b = 4.0 ± 0.3 GeV~ 2 . 
The histogram shows the p 2 t distribution of the DIFFVM Monte Carlo, b) do lp /dp 2 t for J '/if; 
meson production with proton dissociation integrated over 30 GeV < W 7P < 150 GeV. The 
straight line is an exponential fit to the data in the range p 2 < 1 GeV 2 which yields a value of 
b — 1.6 ± 0.3 GeV~ 2 , the dashed line represents the DIFFVM Monte Carlo calculation. 

4.2 J/ip Production with Proton Dissociation 

Energy dependence and p\ distribution 

For the measurement of the proton dissociation cross section the two track sample with a signal 
in one of the forward detectors is used as defined in section 01. The Monte Carlo generator 



DIFFVM for proton dissociation with parameters as in section |3.4| is used for acceptance and 
efficiency determination. The procedure for calculating the cross section and systematic errors 
is the same as in the elastic case. The main contribution to the systematic error is due to 
the selection of the event sample with the help of the forward detectors. The size of the error 
is essentially the same (12%) since both cross sections, the elastic and the proton dissociation 
cross section, are nearly of the same magnitude. The systematic error of the acceptance may be 
slightly larger than in the elastic case due uncertainties in the Mx dependence, this is however 
not taken into account. 

The 7£> cross section for J/if> production with proton dissociation derived from J /if; decays 
to and e + e~ as a function of W 1P is given in table [| and shown in Fig. |3|b. It can be 

fitted by a W* dependence with S = 1.2 ± 0.2 (statistical and systematic error), slightly larger 
than the value for the elastic data, which may be explained by the increased phase space for 
producing excited nucleon systems at higher energy. 

The cross section da ep /dp 2 for J/ip production with proton dissociation is shown in Fig. [| b. 
A log-likelihood fit below 1 GeV 2 gives: 



1.6 ±0.3 ±0.1) GeV 



-2 
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Figure 5: da/dcosd* for diffractive J/if> production (elastic and proton dissociation). The 
inner error bars are statistical, the outer ones contain statistical and systematic errors added 
in quadrature. The curve is a fit of the form oc (1 + const cos 2 9*). 



The first error is statistical while the second one represents the systematic error estimated as 
in the elastic case. The error in the slope parameter due to using p 2 instead of t is calculated 
in the simulation to be —6.5%. The slope parameter is roughly a factor 2 smaller than for pure 
elastic scattering. 



Decay angular distribution 

The polarisation of the J/if> can be accessed via the angular distribution of the J/if> decay 
leptons. The angle 9* is used, which is the angle in the J/ if) rest frame, between the direction of 



the positively charged decay lepton and the J /if) direction in the 7p cms (helicity frame) ||, [43 
Assuming s-channel helicity conservation the expected angular distribution is: 

da 2 1 + cos 2 9* 
oc (1 — p) sin 9 + p 



dcos 9* 



where p, the fraction of transversely polarised J/ip mesons, is predicted to be one. 

The ep cross section for both decay channels J /if) — > p + p~ and J /if) — > e + e~ as a function of 
cos 9* is shown in Fig. |5|. The full two track sample is used in this analysis including elastic and 
proton dissociation samples. A y 2 fit yields p — 1.2 ± 0.2, which is consistent with s-channel 
helicity conservation. 
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Figure 6: Mass distribution for u. + u~ of the inelastic event sample. The curve is a fit of a 
Gaussian plus a polynomial background to the J/ip mass region. The shaded histogram shows 
the mass distribution of like sign muon pairs. The maximum of the fit is at (3.10 ± 0.01) GeV 
with a width of 65 MeV. 



5 Inelastic J/i/j Production 



The event selection for inelastic J / ip production starts from the preselected lepton pair sample 
described in section [3.5| . The mass distribution for fi + fi~ pairs with at least one additional track 
from the interaction point is shown in Fig. |||. The inelastic J/ip candidates are selected by a 
cut around the nominal mass of ±225 MeV. The background is due to muons from leptonic tt 
and K decays and misidentified hadrons. It is subtracted by using the data above and below 
the mass peak as an estimate of the background. 

The elasticity z - calculated according to eqn. || in section || - is used to define the inelastic 
sample. In Fig. |7| an overview over all candidates for J/ijj — > yU + /i~ is shown as function of 
z. The measured z distribution is shown for the inelastic J/ip candidates with z > 0.3. The 
total two track sample and also the subsample with activity in the forward detectors which 
were used in the previous section to evaluate the elastic and proton dissociation cross sections, 
respectively, are shown in the highest z-bin of the figure. The data in Fig. ^ are compared with 
three Monte Carlo models: the contribution of photon gluon fusion as modelled in the generator 
EPJPSI (LO colour singlet model) is shown, at high z the contribution of proton dissociation 
processes as modelled in DIFFVM can be seen, and at low z the contribution of the hadronic 
component of the photon, the "resolved photon process [44}|", is shown as modelled by EPJPSI. 



Note that the simulations are normalized to the data. DIFFVM with proton dissociation is 
normalised to the two track data with forward detector signal. The contribution of photon 
gluon fusion is normalised to the inelastic data for 0.45 < z < 0.90 and the resolved component 
is multiplied by the same normalisation factor. 
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Figure 7: Distribution of the elasticity z for J ftp — > . The HI data in tie J/ ?/> mass 

window are corrected for background below the mass peak. They are classified as events with 
exactly two tracks irrespective of any forward tag (full square), events with exactly two tracks 
and a forward tag (open square) and inelastic events with z > 0.3 (full circle). For comparison 
the following model calculations are shown: a diffractive model including proton dissociation 
(DIFFVM) normalized to the data (two tracks with forward tag); a LO photon gluon fusion 
model (LO BGF EPJPSI) for inelastic J j ip production using the gluon distribution from the 
MRSD- ' parametrisation, and the resolved contribution calculated in EPJPSI with the GRV- 
LO photon structure function. The resolved contribution is scaled with the same normalisation 
factor as the contribution from photon gluon fusion. 



An inelastic sample is selected by requiring 0.45 < z < 0.90 which rejects diffractive and 
resolved photon events and yields a total of ~ 100 J/ip candidates. The data sample is cor- 
rected for efficiency and acceptance using the Monte Carlo generator EPJPSI f29j with the 



parametrisation of the gluon density of MRSD-' [[EJ. The cross section for inelastic J J if) pro- 
duction is determined for < z < 0.90 in order to compare with theory. Therefore a correction 
is applied for the loss below z < 0.45, which is estimated to 13.5 ± 2%, where the error is 
due to using a different gluon distribution (MRSD0'). An additional correction takes into 
account the contamination by J ftp meson production via the resolved photon process. This 
background is estimated using an option of the Monte Carlo generator EPJPSI. Using the 
GRV-LO-parametrisation [46] of the photon structure function yields a background estimate of 
2.5% for z > 0.45. The contribution changes by +2% using as photon structure function the 



parametrisation of LAC1 p7 



Systematic errors for trigger and selection efficiency are estimated as in the elastic analysis. In 
addition to the relevant errors quoted in table |3] which yield a contribution of 12%, a contribution 
to the systematic error of 10% is estimated specifically for the inelastic cross section. This 
number contains the uncertainty due to background subtraction, an error of the contamination 
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J /if) — > 


W w [GeV] 


30-60 


60-90 


90-120 


120-150 


events 

(bg. corrected) 


8.0±3.7 


27.0±6.1 


23.0±8.3 


27.0±7.8 


^-analysis 
^trigger 


0.326±0.013 
0.333±0.027 
0.304±0.024 


0.799±0.032 
0.387±0.031 
0.290±0.023 


0.854±0.034 
0.329 ±0.026 
0.413 ±0.033 


0.731±0.029 
0.279±0.022 
0.646±0.052 


a(ep — > eJ/ipX) [nb] 


1.5±0.7±0.2 


1.8±0.4±0.3 


1.2±0.4±0.2 


1.3±0.4±0.2 


$ 7 /e 


0.0736 


0.0370 


0.0229 


0.0149 


a( 7 p^ J/^X)[nb] 


20.0±9.4±3.2 


49.5±11.1 ±7.9 


52.6±19.0±8.4 


84.1±24.3±13.4 



Table 4: Inelastic J /if) production: Efficiencies and cross sections for inelastic J /if} — > 
production with < z < 0.90. The number of events is background corrected. $ 7 / e is the 
photonflux integrated over Q 2 and W 1P . The errors of acceptance, analysis and trigger efficiency 
are statistical and systematic errors added in quadrature. The hrst error of the cross sections 
is statistical, the second is systematic. 

with resolved events, the correction for the low z region. A further contribution to the error of 
5% is due to a systematic shift in the z determination after reconstruction. In the simulation 
z is found to be systematically shifted to higher values after reconstruction with respect to the 
generated value by 4% for z = 0.8 and by 10% for z = 0.5. 



Energy dependence of ^p cross section 

The 7p cross section for inelastic J /if) production can be derived in a similar way as in section 



4.1. Efficiencies, event numbers and the final cross section for z < 0.9 are given in table 0. The 



resulting 7p cross section is shown in Fig. |8]a. The data are compared to calculations in the 



colour singlet model including NLO contributions |L5 . 

As input the theoretical calculations use the charm mass m c = 1.4 GeV, A-^r = 300 MeV, the 
renormalisation scale is chosen to be identical to the factorisation scale, \f2m c . Predictions are 
shown for different gluon density distributions that have been derived from recent data. The 
gluon distributions differ in the low x behaviour where x is the fraction of the proton momentum 
carried by the gluon. The distributions used can at low x be parametrised as x g(x) oc x _A , and 
the values for A range between (MRSD0') gjj) and 0.4 (MRSG |8|). In Fig. §a the curves 



contain an additional contribution of inelastic if)' production (« 15%, estimated in |15]|). They 
reproduce approximately the energy behaviour of the data, but cover a wide range in absolute 
normalisation. The agreement with the data becomes better with increasing steepness of the 
gluon density at low x. 



The NLO calculation is not fully under control for p\ — > and z — > 1 [48]. Missing contributions 



of even higher order cause problems as can be seen in Fig. |9| where the theoretical p\ distribution 
bends over for p\ £ 1 GeV 2 . 

Therefore the data and prediction are shown in Fig. [8|b for a restricted kinematical range, 
z < 0.8 and p\ > 1 GeV 2 and thus requiring the emitted gluon to be hard (thereby reducing the 
data sample by a factor 2). Improved agreement is observed with all gluon density distributions 
but the sensitivity is reduced. 
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Figure 8: Total cross section for inelastic J / ip -production for z < 0.9 (a) and z < 0.8 andpj > 
1 GeV 2 (b). The inner error bars of the HI points are statistical, the outer include statistical 
and systematic errors added in quadrature. The curves represent NLO calculations for 
different gluon distributions and contain a 15% correction taking into account ip' background. 
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Figure 9: a) da lp /dp 2 for inelastic J/ip production with z < 0.9 integrated over 30 GeV < 
W 1P < 150 GeV. The straight line is an exponential fit to the data in the range p 2 t < 5 GeV 2 , 
the dashed line is a NLO calculation Jl5|/ at W 1P = 100 GeV with the GRV structure function, 
b) da rp /dz for inelastic J/ip production with a cut pf > I GeV 2 , in comparison with NLO 
calculations in the colour singlet model with the GRV structure function. The LO colour octet 
calculation is from 
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p 2 , 2-distribution 

In Fig. |9|a the distribution of the transverse momentum p\ of the J/ip is shown for z < 0.9. A 
log-likelihood fit to the data of an exponential exp (—bp 2 ) in the region p\ < 5 GeV 2 yields: 

b = (0.39 ± 0.06 ± 0.03) GeV" 2 
where the first error is statistical and where the systematical error is estimated by using a a 



X -fit and varying the fit region. This result is in agreement with the NLO-calculation [|15 
which predicts a slope of b = 0.3 GeV~ 2 above a p\ > 1 GeV 2 . 

The differential energy distribution of the J/ip, da/dz, is shown in Fig.|9|b with a cut in p\ > 
1 GeV 2 . It is compared with the NLO calculation (ff5j in the colour singlet model with the 



GRV [|39| structure function for the proton. Agreement in shape and normalisation is found 
within errors. This is interesting in view of the speculations about possible additional colour 
octet contributions to photoproduction of J/ip. These yield a large inelastic contribution at 
z > 0.8 |17] after fixing the normalisation by the CDF data|l6|. The strong increase indicated 
in Fig. R|b is not supported by the present data. 



6 ifj' Production 

The same data sample as for the J/ip analysis is used to measure the cross section for the 
diffractive photoproduction of ip' mesons (i/j(2S)), by searching for the decay ip' — > J/ipT[ + 7T~, 
where the J/ ip subsequently decays to a muon pair. 
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For this investigation, only the muon triggers are utilised, but the identification of the muon 
pair is otherwise identical to that used in the J/ip analysis. It is further required that, in 
addition to the muon pair, exactly two other, oppositely charged tracks are found in the central 
tracking detector which are assumed to be pions. For events in which both pion candidates 
are measured in the central tracking detector to have a transverse momentum greater than 150 
MeV, a mass difference is formed, defined as AM = ra lJL + ll -^+ n - — rn^+^-. A peak of 7 events 
is found in an interval of ±60 MeV around the mass difference, AM = m^> — mj/^, with no 
background events seen outside the peak. 

Monte Carlo studies show that this selection restricts the kinematic domain of acceptance to 
z' > 0.95, 40GeV < W 1P < 160 GeV and Q 2 < 4GeV 2 , where z' = E^/E^ in the proton rest 
frame. The acceptance is estimated using the same Monte Carlo generator as was used for 
the J/ip analysis, with parameters of the model chosen to match the J/ip data. A mixture of 
elastic and proton dissociation events is assumed, which is compatible with that found in the 
J/ip events. No separation of elastic events from events with proton dissociation is attempted. 
In the kinematic region above, the combined acceptance of trigger and selection is found to be 
4.4±0.9%. Assuming a branching ratio for this ip' decay of 0.019±0.002 p[9" |, the ep cross section 
for the production of ip' in the kinematic domain z' > 0.95, 40 GeV < W yp < 160 GeV and 
Q 2 < 4 GeV 2 is 2.9 ±1.1 ±0.6 nb, where the first error given is statistical, the second systematic. 
This translates into a photoproduction cross section, cr 7P (W 7P = 80 GeV) = 24 ± 9 ± 5nb for 
ip' production with z' > 0.95 and Q 2 < 4 GeV 2 . The corresponding cross section for J/ip 
production is the diffractive cross section (including elastic and proton dissociation) which at 
W yp = 80 GeV is a 7P = 119.5 ± 11.2 ± 12.0 nb. The ratio of ip' to J/ip production is then 
0.20 ± 0.09 which is in agreement with previous measurements |9|, |4l], |5( 



7 Summary 

J/ip meson production in the photoproduction limit is analysed in the elastic channel. Correc- 
tions for contributions from events with proton dissociation are applied almost entirely on an 
event to event basis. The following results are found: 

• The total jp cross section is observed to increase with energy as W^ with 5 = 0.64 ±0.13 
for HI data alone, and 5 = 0.90 ± 0.06 including the ZEUS and low energy data. This 
increase is faster than predicted in soft diffractive models (5 = 0.32). 

• A calculation using the model by Ryskin in the framework of perturbative QCD results in 
a good description of the energy dependence if the gluon distribution from the set MRSA' 
is used. 

• The p 2 distribution of the elastic J / ip data below 1 GeV 2 can be fitted with an exponential 
e"*?, b = 4.0 ± 0.2 ± 0.2 GeV -2 . 

The 7p cross cross section for J/ip meson production with proton dissociation was determined. 
The 7p cross section is as large as the elastic one and its energy dependence is slightly steeper 
than the elastic cross section, 5 — 1.2 ± 0.2. The p 2 distribution is flatter than for the elastic 
process. The slope parameter is approximately a factor 2 smaller than for elastic data. 
The angular distribution in the helicity frame was determined for the complete diffractive data 
sample, i.e. including elastic and proton dissociation contributions. The distribution in the 
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helicity frame is consistent with s-channel helicity conservation, p — 1.2 ± 0.2, where p is the 
fraction of transversely polarised J/ip mesons. 

Inelastic J/ip meson production can be described well by a QCD calculation in NLO of the 
colour singlet model if cuts are used which demand the emitted gluon to be hard, pt > 1 GeV 
and z < 0.8, in order to make perturbation theory applicable. The agreement is observed in the 
shape of the W 1P , pf and z distributions. The absolute normalization between data and theory 
agrees within the experimental errors of m 30%, using any gluon distribution function which 
describes the F 2 data at low x. A comparison to the relative energy distribution z of inelastic 
J /if) production with pf > 1 GeV 2 calculated in a colour octet model shows disagreement at 
high values and makes a large colour octet contribution unlikely. 

Diffractive production of ip' was observed. An estimate of the ratio of the cross section to the 
one for J/tp yields 0.20 ± 0.09 in agreement with previous measurements at lower energy. 
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